Dusty Mgll Absorbers: Implications for the GRB/Quasar 

Incidence Discrepancy 



Vladimir Sudilovsky^'^ 
Donald Smith^ 
Sandra Savaglio^ 

ABSTRACT 

There is nearly a factor of four difference in the number density of intervening 
Mgll absorbers as determined from gamma-ray burst (GRB) and quasar lines of 
sight. We use a Monte- Carlo simulation to test if a dust extinction bias can 
account for this discrepancy. We apply an empirically determined relationship 
between dust column density and Mgll rest equivalent width to simulated quasar 
sight-lines and model the underlying number of quasars that must be present to 
explain the published magnitude distribution of SDSS quasars. We find that an 
input Mgll number density dn/dz of 0.273 ± 0.002 over the range 0.4 < z < 
2.0 and with Mgll equivalent width Wq > I.OA accurately reproduces observed 
distributions. From this value, we conclude that a dust obstruction bias cannot 
be the sole cause of the observed discrepancy between GRB and quasar sight- 
lines: this bias is likely to reduce the discrepancy only by ~ 10%. 

Subject headings: surveys, dust extinction, quasars: general 



1. Introduction 

The Sloan Digital Sky Survey (SDSS) has greatly expanded our knowledge of the large 
scale structure of the universe, in large part by greatly increasing the nunijper o f known 



quasars at high redshift to ~ 10 (iPetersonI 119971 : lAdelman-McCarthy et al.ll2008l ). These 



bright, distant sources of light also provide a means to probe intergalactic space by al- 
lowing us to identify and analyze spectral absorption lines caused by intervening matter. 
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These measurement s can constrain the star formation rate history (jHamann fc Ferlandlll999 



Pettini 



Hopkins et ah the formation of large scale structure (IClowes et al. I ll999l:l Cabanac et al 

2005!) . the evolution of chemical enrichment in the high -redshift universe teaker et al lboool: 



2001 



Mehlert et al 



20061 : iMesinger fc HaimanlboOTh . 



2003 : Prochaska et al. and the time of re-ionization (iFan et al.l 



These investigations are only made possible by a population of bright sources at high 
redshifts. Quasars provide such a population, but so do the brief flares of Gamma-Ray 
Burst (GRB) sources. The catalog of quasar observations is much larger than the sample of 
GRBs that have been studied with optical spectroscopy, in large part due to the transient 
nature of GRBs. GRBs are first detected in gamma-rays, thus they are therefore not subject 
to the same observational biases as quasars. For this reason, the GRB population plays a 
potentially crucial role in the study of intervening absorption systems. 



Prochter et al.l (120061 ) used GRBs to study absorption systems that lie along the hues 
of sight. They found on average four times the number of intervening strong (equivalent 
width greater than one angstrom) Mgll absorbers in GRB sight-lines per unit redshift as 
compared to the number of Mgll absorbers in quasar sight-lines. This result is surprising 
because both populations sample random, independent sight-lines. Possible factors that may 
explain this discrepancy include differi ng Mgll covering factors in GRB and quasar beams 
( iPontzen et al.l 120071 : iFrank et al.ll2007l ). weak gravitational lensing affecting a population of 
absorbers, a dust extinction bi as, and host-associated Mgll absorbers bei ng incorrectly iden- 



tified as intervening systems (jProchter et al.l l2006l : iPorciani et al.l l2007l ) . ISudilovsky et al. 



( I2OO7I ) surveyed intervening CIV systems and found no difference in th ei r number density 



i n eith er type of sight-line, which is also confirmed by iTejos et al.l (120071 ). ISudilovsky et al. 



( I2OO7I ) argued that the simplest explanation for this discrepancy hes in a selection bias. 
Samples of quasars that are detected based on optical magnitudes may be biased towards 
sight-lines that contain fewer intervening Mgll systems, because Mgll systems are tracers of 
dense, potentially dusty c l ouds, which wou l d attenuate optical hght passing through them 
dMenard et aPboOsl . boosi : IVork et al.l[2006l : iKhare et al.1l2005bh . 



We report on the results of a Monte-Carlo simulation to test the selection bias introduced 
by dust extinction on optically limited surveys. We simulated a population of quasars that 
follow distributions in redshift and magnitude determined empirically from the SDSS. We 
applied a magnitude shift to each quasar that corresponds to the effect of the dust associated 
with simulated Mgll absorbers placed along that quasar's line of sight. For various input 
Mgll number densities dn/rfjl], we added quasars to the simulation until the number of 



For the remainder of this paper, we constrain the quantity dn/dz to 0.4 < z < 2.0 and Mgll equivalent 



3 



quasars brighter than the SDSS's hmiting magnitude matched the number of quasars in 
the actual SDSS catalog. Finally, we determine the input dn/dz for which the simulation 
returns the SDSS-observed distributions of Mgll dn/dz and quasar i-magnitude. Using this 
technique, we determine if dust can explain the GRB and quasar Mgll number density 
discrepancy. 

In this paper, we explain the observationally-derived probability distributions upon 
which our simulation is based in §21 In §3] we present the algorithm used by our simulation, 
and in §1] we present the statistical properties extracted from many simulation runs. §3 
places these results in the context of current research. 



2. Observational Data 



2.1. The Mgll Sample 



Mgll absorption is readily detecte d in the optical spectra of objects with 0.4 z < 
2.2 (ILanzetta. Wolfe, fc Turnshekl 119871 ). Mgll is fo und at small galaxy impa c t parameters 



and can therefore be used as a tracer of galaxies (ILanzetta fc BowenI Il990l : ISteidel et al. 



19941 ). The optical depth of dust has been well approximated over 0.4 < z < 2.0 and Mgll 
equivalent width Wq > l.OA with a power law of the form (xy)(l^o) = '^v,o (■^l , where 



a 



1.8 8 ±0.17, r„,n = ('2.5 ± 0.2) x 10"^ and Wo is the rest equivalent width of the A2796 



feature (IMenard et al.ll2008l ). 



We calculated the extinction caused by dust associated with Mgll absorbers as a function 
of their equivalent width Wq, their redshift, and t heir a bundance along sight-lines to distant 
objects such as quasars or GRBs. iMenard et al.l (120081 ) determined the correlation between 
Mgll Wo and the optical depth of dust, and the dust optical depth yields the extinction 
calculation for one absorber system. We limited our simulated absorbers to Wo > 1 A 
to remain consistent with the domain of applicability of the Mgl l number density and Wp 
distrib utioij^. We use the Wo distribution empirically derived by Nestor. Turnshek. fc Rao 
( I2OO5I ). Since it is possible that there exist a population of strong Mgll absorbers that 
contain so much dust that they simply cannot be observed through optical spectroscopy, our 
estimate of the fraction of dust-obscured quasars should be regarded as a lower limit. We use 



width Wo > l.OA 

^Also, evidence suggests that dust extinction is insignificant below this cutoff ( York et al. 20061 : 



Menard et aUboOSh 
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the Mgll dn/dz as determined by lProchter et al.l (120061 ) as an upper limit in the simulation 



2.2. The Quasar Sample 

We used the main spectroscopic quasar sample from the SDSS fifth data release (DR5) 
as our control catalog. The DR5 differentiate s quasars from st ars by color and redshift 
criteria according to an algorithm described by I York et al.l (120001 ) . The SDSS main quasar 
sample consists of quasars in the magnitude range 15.0 < i < 19.1. We exclude quasars with 
i > 19.1 from our analysis. These quasars are detected by the SDSS serendipity fibers and 
are therefore not detected using a well-defined search algorithm. 

We further constrained this quasar catalog to the redshift range z < 2.0. As described 
in §2.H the number density of intervening Mgll systems is well-constrained only within this 
redshift range. This left us with a final sample size of Ngdss = 11361 quasars. The i-band 
and redshift distributions of the SDSS sample can be found in Figure [TJ 



3. Analysis Methods 

We produced a sample of simulated quasars that has the same statistical properties as 
the original SDSS sample. For every simulated quasar, we added a number of Mgll absorbers 
whose characteristics are based on observationally-determined distributions to the sight-line. 
We accounted for the effect of d ust by dimming the i -band magnitude of the quasar according 
to a SMC-type extinction law (jPrevot et al.lll984j : iMenard et al.ll2008l ). Evidence suggests 
that SMC-like dust is the most comnaon form of dust associated with intervening Mgll 
dMenard et aPboOsI : iKhare et al.lboOSal : Iwild et al.l[2006[ ). 



We applied 25 different input values for Mgll dn/dz in the range 0.24 < {dn/dz)in < 
0.90. For a given {dn/dz)in, we added quasars to the simulation until the number of quasars 
whose i-band magnitude was brighter than 19.1 equaled the number of quasars in the SDSS 
sample, l^sdss- Once this stopping criterion was met, we were left with a sample of N^j^ 
simulated quasars. Since extinction dimmed some of these quasars below 19.1, N^j^ was 
always larger than l^sdss- We interpret the quantity sim—'^ sdss as the number of quasars 
the SDSS misses due to dust dimming the i magnitudes past the i=19.1 limit. We repeated 
this step 500 times and recorded the mean number of quasars in the simulation (Nsim), its 
associated uncertainty cjv^.^, the mean output Mgll number density {{dn/dz)out) , and its 
uncertainty (J(dn/dz)out^ all for the given [dn/dz)in. Finally, we quantified the relation between 
intrinsic and observed Mgll number density and its associated uncertainty. 
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4. Results 

For a given {dn/dz)in, we created quasars until the stopping criterion described in ^ 
was met. The primary result from a single run of the simulation was a sample of quasar 
magnitudes. Several runs' samples for various values of {dn/dz)in are plotted in Figure [2] 
as a solid line histogram, along with the i-band distribution of quasars in the main SDSS 
catalog as a dashed line histogram. The same total number of quasars are brighter than 
19.1 in each histogram. The solid bins to the right of i=19.1 represent the simulated quasars 
missed by the SDSS due to dust effects. 

From a single run with a given {dn/dz)in, we derive the value for the Mgll number 
density that a hypothetical observer would deduce from observing the populations of quasars 
that survive the i < 19.1 cut {{dn/dz)out)- We show {{dn/dz)out) as a function of {dn/dz)in in 
Figure[3l where the inset highlights the data in the {dn/dz)in ~ 0.25 range. By interpolating 
the data with a cubic polynomial, we find that {{dn/dz)out) = 0.24 when {dn/dz)in = 
0.273 ± 0.002. We estimate the uncertainty in this value by re-running the simulation with 
values of {dn/dz)in between 0.265 and 0.281 and finding the values for which {{dn/dz)out) 
were different from 0.24 95% of the time. 

We present the distribution of Ngi^ — Ngdss for an input number density of 0.273 in 
Figure m The data follow a Gaussian distribution, with (A^sim — N^dss) = 280 ± 20 quasars. 
This implies that the SDSS misses (2.4 ± 0.2)% of 2; < 2.0 quasars due to dust extinction 
associated with intervening Mgll absorbers. 



5. Discussion 



According to the simulation, the discrepancy between Mgll number densities seen in 
GRB and quasar sight-lines cannot be fully explained by a selection bias introduced by dust 
extinction. We find that the observed incidence of intervening Mgll absorbers in quasar 
sight-lines is lowered due to dust extinction by ~ 10%: The unbiased Mgll dn/dz is likely 
to be 0.273, which is still very different from the GRB-observed value of dn/dz = 0.90. 
Additionally, we find that at least ^ 2% of z < 2.0 quasars are not included in the SDSS 
main quasar s ample due to dust effects. This value is consistent with the estimate by 
Menard et al.l (120081 ) in which they use a different approach. 



Porciani et al.l (120071 ) have also tested explanations for the GRB and quasar Mgll num- 



ber density discrepancy. Most relevantly to t his work, they use mock light-cones drawn 
from the Millennium Run (jSpringel et al.ll2005l ) to model the number of galaxies that would 
obscure a population of quasars at z = 2.3, assigning Mgll absorbers probabilistically to 
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galaxies and assuming a relation between Mgll and dust that best reproduces observed col- 
ors. Using this technique, they conclude that 16% of quasars are missed by the SDSS due 
to dust effects, and that dust extinction can account for the GRB/quasar number density 
discrepancy only by a factor of ~ 1.3 — 2. A significant difference between this work and 
the work done by Porciani et al. is in how we calculate the dust traced by Mgll. Porciani 
et al. use a Weibull distribution for color excess E(B-V) values, and assign these values to 
intervening Mgll absorbers. A Weibull distribution is a continuous probability distribution 
whose shape is determined by two parameters. Besides constraining these two parameters 
with the quasar relative color excess A{g — i), the E(B-V) distribution is a free parameter 
in their model. In this work, we use an empirically derived relationship between Mgll and 
dust reddening. Additionally, our method does not rely on the semi-analytical modeling in 
the Millennium Run. 

Though a bias stemming from dust extinction is undoubtedly present in the SDSS, it 
only accounts for a small fraction of the observed Mgll number density discrepancy between 
GRB and quasar sig ht-lines. Furthe r more , none of the other proposed solutions to this prob- 



lem are satisfactory: iPontzen et al.l (120071 ) show that there is no evidence for systematically 
smaller Mgll equivalent widths over quasar broad emission regions (as compared to quasar 
continuum regions) and conclude that di fferent Mgll cove r ing fa ctors in the two types of lines 



of sight cannot explain the discrepancy. ICucchiara et al.l (120081 ) investigate the properties of 



Mgll absorbers in GRB and quasar lines of sight and find no difference, there by arguing that 



the ab sorbers are correctly identified as intervening, non-intrinsic systems. iPorciani et al. 



( 120071 ) point out that a bias stemming from gravitational lensing requires that the quasar 
beam be larger than the GRB beam and that the optical depth of micro-lenses is likely 
significantly greater than observed. Both of these conditions, though possible, are unlikely 
given current observational data. 

Although the Mgll number density discrepancy is statistically significant, it is neverthe- 



l ess im portant to increase the GRB absorber sample size to better constrain it. IChen et al. 



( 120091 ) have found galaxies within 2"of GRB hosts whose afterglows showed strong Mgll 



absorbers. A comprehensive survey of these galaxies may also shed new light on this issue. 

In summary, we have used a Monte- Carlo simulation to test if the factor of ~ 4 difference 
in Mgll number densities observed in GRB and quasar sight-lines is due to dust extinction. 
We find that a dust extinction bias is likely to account for only ^ 10% of the observed 
overdensity in GRB sight-lines. Additionally, we estimate that ~ 2% of 2; < 2.0 quasars are 
not included in the SDSS main sample because of dust obstruction. 
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Fig. 1. — Left: Redshift distribution of Nsdss quasars in bins of Az — 0.0046. Right: i-band 
magnitude distribution of Nsdss quasars in bins of Ai = 0.0096. The y-axis represents the 
number of quasars in each bin. These quantities have been extracted from the SDSS DR5 
catalog. 
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Fig. 2. — Dashed line: ^-band magnitude distribution of the observed f^sdss quasars. Solid 
line i-band magnitude distribution of simulated N^jm quasars. Example simulation outputs 
are plotted for {dn/dz)in — 0.24,0.273, and 0.90, respectively. 
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Fig. 3. — {{dn/ dz)out) as a function of {dn/dz)in, with the sohd hne representing equahty. 
The data are well fit by a cubic polynomial. The inset highhghts the data in the range 
0.24 < {dn/dz)in < 0.30. {{dn/dz)out) = 0.24 when {dn/dz)in = 0.273. 
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Fig. 4. — Distribution of Ngdss — Nsimi which represent the number of quasars in the sample 
that would be observable if there were no dust extinction associated with Mgll absorbers. 
The distribution is Gaussian, with {Nsim — Nsdss) = 280 ± 20 quasars. 



